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The combined diffusion coefficient method is widely used to treat the mixing and demixing of different 
plasma gases and vapours in thermal plasmas, such as welding arcs and plasma jets. It greatly simplifies the 
treatment of diffusion for many gas mixtures without sacrificing accuracy. Here, three subjects that are 
important in the implementation of the combined diffusion coefficient method are considered. First, it is 
shown that different expressions for the combined diffusion coefficients, arising from different definitions 
for the stoichiometric coefficients that assign the electrons to the two gases, are equivalent. Second, an 
approach is presented for calculating certain partial differential terms in the combined temperature and 
pressure diffusion coefficients that can cause difficulties. Finally, a method for applying the combined 
diffusion coefficients in computational models, which typically require diffusion to be expressed in terms of 
mass fraction gradients, is given. 

he combined diffusion coefficient method 1 allows major simplifications in the modelling of thermal plasmas 
in gas mixtures, in most cases without loss of accuracy. Instead of solving conservation equations for the 
mass of each species (e.g., N 2 , N 2 + , N, N + , N 2+ , N 3+ , Ar, Ar + , Ar 2+ , Ar 3+ , e~ in an argon-nitrogen plasma), 
the species are combined into their parent gases, so only one conservation equation (e.g., for argon or nitrogen 
'gas') has to be solved. The multicomponent diffusion coefficients for each pair of species (of which there are 
q(q — 1)12 ordinary and q — 1 thermal diffusion coefficients for a plasma containing q species) are combined into 
just four combined diffusion coefficients, one each describing diffusion due to composition gradients, temper- 
ature gradients, pressure gradients and applied electric fields. 

The method has been successfully applied to modelling a wide range of phenomena in thermal plasmas. These 
include demixing in welding arcs 2,3 , plasma torches 4 , arc plasma reactors 5,6 , and inductively-coupled rf plasmas in gas 
mixtures 710 , mixing of metal vapour in welding arcs 1113 , mixing of air into plasmas 14 23 , mixing of polymer vapours 
into the arc in circuit breakers 24 29 , and expansion of metal vapour ablated by a laser into the surrounding air 30 . 

Values of the combined diffusion coefficients have been presented for many gas mixtures, including argon- 
nitrogen 31 , argon-oxygen 31 , argon-air 32 , argon-helium 33,34 , argon-hydrogen 35 , nitrogen-hydrogen 36 , sulfur hex- 
afluoride-tetrafluoroethene 29 , argon-metal vapours 11,37,38 , and air-metal vapours 11,3 ' 40 . The methods have also 
been extended to two-temperature plasmas 41 " 43 . 

Although determination of the combined diffusion coefficients is a relatively simple extension of the calcula- 
tion of the multicomponent diffusion coefficients, three points of confusion have frequently been apparent in 
communications with researchers performing such calculations. The first concerns the method of assigning 
electrons to the two gases using stoichiometric coefficients, the second the treatment of the mole fraction gradient 
terms in the calculation of the combined pressure and temperature diffusion coefficients, and the third the 
method of application of the combined diffusion coefficients in computational models, which typically require 
the diffusion term to be written in terms of mass fraction gradients. 

The purpose of this paper is to clarify these three points. To put this in context, an introduction to the combined 
diffusion coefficient approach is first presented. 

In the most general treatment of thermal plasmas containing more than one chemical element, the mass 
conservation equations for each individual species (e.g., Ar, Ar + , Ar 2+ , N 2 , N 2 + , N, N + , N 2+ , e~, etc., in an argon- 
nitrogen plasma): 

^ + V-(pvY i ) + V-J i = r i , (1) 

have to be solved everywhere in the plasma. Here p is the mass density, v is the mass-average velocity, T, is the 
mass fraction of species i, and r, is the net rate of production of species i due to chemical reactions, vaporization, 
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etc. For a gas or plasma containing q species, the diffusion mass flux 
of species i, is given by 44 ' 45 

2 q 

/, = rriitiiVi = 22 mjDijdj - DjV In T, (2) 

P ;=1 

where v, is the diffusion velocity of species i (relative to the mass- 
average velocity), m,- and n { are respectively the mass and the number 
density of the species i, and d, describes the driving forces for dif- 
fusion of species j. The ordinary diffusion coefficients D it and thermal 
diffusion coefficients Dj are in this case multicomponent diffusion 
coefficients 44,45 . 

For thermal plasmas, in which values of El n (electric field normal- 
ized to the gas number density) are very low (typically less than 1 Td 
= 1 X 1CT 21 Vm 2 ) and the species velocity distributions can be 
closely approximated by a Maxwellian distribution, it is standard 
practice to use the Chapman-Enskog approximate solution to the 
Boltzmann equation to obtain the diffusion coefficients. Their cal- 
culation requires values of the mole fractions and the masses of all the 
species present and the collision cross-sections for binary interac- 
tions between each pair of species present 45 . Note that the much 
higher values of Eln present in swarms and many non-equilibrium 
discharges require different approaches 46,47 . 

The driving force term dj is given by 44,45 

1 



VlnP-^-M — F,— Vn,F ( 
Pp \m, > 



(3) 



The three terms describe respectively diffusion due to gradients in 
the mole fraction x p the pressure P, and the external forces Fj acting 
on species j. Here the only external force that is considered is that due 
to an applied electric field. It is interesting to note that gravity does 
not exert a driving force, since the third term cancels to zero. 

Because diffusion velocities are defined with respect to the mass- 
average velocity, only N — 1 need be calculated, and (1) has to be 
solved for only N — 1 species. In plasmas, the diffusion coefficients 
have to be modified to take into account ambipolar diffusion 1,45 . This 
arises because electrons diffuse more rapidly than ions because of 
their lower mass, inducing an electric field that accelerates the ions 
and slows the electrons. 

Clearly, solving the conservation equations (1) for each species, 
and calculating ordinary diffusion coefficient for each pair of species 
and thermal diffusion coefficients for each species, is computation- 
ally expensive. As a consequence, simplified methods are generally 
used; a range of methods of varying levels of simplicity and accuracy 
is available 48,49 . The combined diffusion coefficient approach has the 
advantage of being mathematically equivalent to the full multicom- 
ponent treatment if three conditions are satisfied: the gases are 
homonuclear, they do not react with each other, and the plasma is 
in local chemical equilibrium (LCE). Although different species will 
usually have different diffusion velocities, the imposition of LCE 
means that the local composition depends only on the local temper- 
ature, pressure and the relative concentrations of the two gases. For 
example, although nitrogen molecules and atoms diffuse at different 
rates, dissociation and recombination reactions will maintain the 
composition of nitrogen gas in accordance with LCE. 

The diffusion mass flux of gas A is 

2 

Ta = - mXm (p^VxE + WbV ln p +Wb e )-WbV ln T > ( 4 ) 

where mX and mg are the average masses of the heavy species of the 
respective gases 

P j P 

mX=12 mkXk / E*fr 

k = 2 / k = 2 

. , , ( 5 ) 



m B - 



-- ^2 m k x k 

k=p+l I 



12 Xk 

k=p+l 



and x B is the sum of the mole fractions of the species of the gas B: 

P 

xj= y~](i+z k )x k 

k = 2 

(6) 

XB= 12 ( 1 + Z k) x k = l-XA' 

k=p+l 

where Z k is the charge number of species k, and is used to divide the 
electrons between the two gases. 

The combined ordinary diffusion coefficient D\ B , combined pres- 
sure diffusion coefficient D^ B , combined electric field diffusion coef- 
ficient -D^ G and combined temperature diffusion coefficient 
describe, respectively, diffusion due to mole fraction gradients, gra- 
dients of the total pressure, externally-applied electric fields and 
temperature gradients. They are calculated by writing Ja and Jg in 
terms of the total mass flux of the species of gas A and gas B respect- 
ively. Electrons arising from ionization of gas A species are assigned 
to gas A, and likewise for gas B, using the stoichiometric coefficients 
kj and k,: 



(7) 



7 B = fci/l+ J2 k Ji~- 
i=p+l 



i=l,i=p+l 



where 



fa =12 ZkXk / xu 



k = 2 
9 



fa = 12 ZkXk / *!' 
k=p+l I 



(8) 



The terms fc, and fc, are used to apportion the electrons between gas A 
and gas B according to the degree of ionization of the respective gas. 
Using 

Vx, = {dx,/dxB)VxB+ (dxj/dP)VP+ (dXj/8T)VT, (9) 

which is valid since our assumption of LCE means that an equilib- 
rium composition can be defined as a function of temperature, pres- 
sure, and the relative concentrations of the two gases, we obtain 1,50 : 



1 ox 



and 



D T AB - 



I P i 
— Y^km^mjXjZjDy, 



7 = 1 



OX, 



(10) 



(11) 



(12) 



(13) 



We can obtain analogous expressions for D^ A , P^ BA , D^a an ^ D B a mat 
satisfy D*J, = - D^, L% = -D^, and ti^ = - D^. 
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Results 

Stoichiometric coefficients. In previous publications 1,3,32,35 ' 48,49 , 
different expressions for the combined diffusion coefficients have 
been given. The derivation of the coefficients was performed using 
number fluxes rather than mass fluxes, leading to the following 
expressions: 

I P <f 
,= — Vs, Vm,D° 



m B - 



dxi 



(14) 



which can be shown to be equivalent to (10) using equations (5) and 
(8). 

Calculation of partial differentials in the expressions for combined 
diffusion coefficients. Equation (9) is more formally written as 

Vx j = (8x j /m) \ P T m+ (8x } /dP) |^ r VP+ (Bxj/dT) |^VT,(21) 

denoting that, for example, [pXj/ '3T), has to be calculated with xg 
and P held constant. Hence, equation (13) is written 



,=1 ;=1 



fa] 

dp 



j p q 



and 



n Ta r? J-^ 



8T 



(15) 



(16) 



(17) 



(18) 



where the s, are stoichiometric coefficients, defined by 



si = — > , 

Wl A ^ 



s t = b^ Xk I ^ b k x k , 2<i<p, 

k = 2 I k = 2 
1 I 1 

Sj = bj ^ x k I £ b k x k ,p+l<i<q 

k=p+\ I k=p+l 



where b f is the number of atoms in a molecule of species i; thus b t = 2 for 
a diatomic molecule or ion, and fc, = 1 for an atom or monatomic ion. 

In the initial publications, Si = 0 was used 1,48 . This is equivalent to 
neglecting the mass flux of the electrons Ji = m\tl\V\, which is 
reasonable when calculating D^ B , D^ B and D\ B , since the electron 
mass is much smaller than the heavy species mass, and the electron 
diffusion velocity is of the same order of that of the heavy species due 
to ambipolar effects. However, ambipolar effects do not play a role in 
electric field diffusion, meaning that V\ » Vj(i j= 1), so Ji cannot be 

neglected in calculating D^g 3 - 

The expressions for the combined diffusion coefficients given in 
equations (14) to (17) are in fact equivalent to the current expressions 
(10) to (13). We show here that the expressions given in equations 
(10) and (14) for are equal; the expressions for the other coeffi- 
cients can be shown to be equal in a similar manner. 

Substituting equation (18) into equation (14) gives 



nil i^ZiXiir 1 -^ OXj 1 , k = 2 a ox i 



D T AB - 



8x- 



(22) 



and similarly for equation (11). 

For a given gas mixture, however, xg is a function of temperature 
and pressure, due to the occurrence of dissociation and ionization 
reactions. In standard computer programs for the calculation of 
plasma composition, the relative mole fractions of the gases are set 
at a given temperature and pressure, typically 298 K and 1 atm. 
Maintaining these mole fractions constant does not ensure that x~b 
stays constant with changes in temperature and pressure. 

This problem can be solved by obtaining expressions for 
(dxj/8T^\— p and (pXj/ 8P)\ — T in terms of the value of xb at a 
particular temperature and pressure. These can be derived using 
the identity 



(23) 





8f 


8x 




dy 


dt 


s 8x 







where / = fix, y), with x = x(s, t) and y = y(s, t). 
If x = r, so y = y(s, x), this gives 



8f 






4 


dy 


ex 


, dx 


y 


dy 


x 8x 



(24) 



Substituting x-. iorf, T for x, x% for y, and x 2 98 k for s » we obtain 

( Sx ,/ BT ) kp = (^/8T) |_ p - (8x,/m) \„ >T x (Sxj/ST) \^ P (25) 

where x 2 9sk is the value of xjs at 298 K. 
Similarly, we obtain 

( 8x j/ dP ) Lt = ( dx i/ 8p ) \^,t ~ (%/^) \p,t x (^/8P)\^t (26) 

where X\ atm is the value of at 1 atm. 

The values of (pXj/dT) in equation (13) and (dXj / oP) in equation 
(11) are then calculated using equations (25) and (26), respectively. 

Figure 1 shows a comparison the combined temperature diffusion 
coefficient calculated with and without the 'correction term', i.e., 
with and without the second term on the right-hand side of equation 
(25). The collision integrals and methods used were those given in 31 , 
as updated in 36 . It is clear that the 'correction term' is very significant. 
The corresponding term in equation (26), used in calculating the 
combined pressure diffusion coefficient, is in contrast, negligible, 
because xjs depends only weakly on pressure. 



Since b t lbj = mj nij for species derived from the same gas, this can be 
written 



mi 



E- ZiX; „ v 

— > m-.If.i- 



f* _j 

] xls mi 



P 

S X k p q 

-jr 1 — 

J2 m kXk 1=2 > =1 

k=2 



?1( 2 °) 



Application of combined diffusion coefficients in computational 
models. The combined diffusion coefficient approach allows us to 
calculate J a, the diffusion mass flux of gas A. The conservation 
equation for gas A, replacing the conservation equations (1) for 
each species i, can be written. 



8pY A 
8t 



+ V-(pvY A )=-V-J A +S A 



(27) 
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where Ya is the sum of mass fractions of the species making up gas A 
and S A is the net rate of production for gas A (for example by 
vaporization of a surface). 

For use in standard computational algorithms, J a has to be written 
in terms of Ya- We can modify Eq. (4) to read 

Ta=-p"^(Wb VxX-WbV^P-Wb *)-DLVlnT (28) 

noting that VxX = — Vxg and n 2 / p = p j M 2 , where M is the average 
mass of all species in the gas mixture 

1 

M=^Tm k x k . ( 29 ) 

k=l 

Since Ya = [Ma / M) *a > where Ma is the average mass of all species of 
gas A: 

P j P 

M A = Y^(mx k + m 1 Z k x k ) r^(\+Z k )x k , (30) 



k = 2 



k = 2 



we have 



VxX=^VY A + Y A V 
M a 



\M A J 



(31) 



Therefore, 



— mAing — mirtin / M 

lA — p im^ YA - p i^^ YA \WA 



+ P — 2 Dab^P + P — D>i„E-Di B V\nT 



(32) 



M 



M 



Substituting this into equation (27) gives 
- + V-(pvY A ) = V- ( p _ _ DabVYa 



ct 



M A M 
m A m B 
M 



M 

Ma 



-V^^Z^lnP 
+ V-fDjIV In T)+S A 



(33) 

rnXmBj^- 
P^=2~DabE 
M 



The first term on the right-hand side is the diffusion term, while the 
second, third, fourth and fifth terms on the right-hand side are not in 
VYa, and therefore have to be treated as source terms. 

Discussion 

The following points related to the combined diffusion coefficients 
have been clarified. 

1. Different expressions for the combined diffusion coefficients, 
arising from different definitions for the stoichiometric coeffi- 
cients that assign the electrons to the two gases, have been shown 
to be equivalent; 

2. An approach to dealing with the partial differential terms for 
which the mole fractions of the gases are held constant, which 
have to be calculated in determining the combined temperature 
and pressure diffusion coefficients, has been presented; 

3. A method for the application of the combined diffusion coeffi- 
cients in computational models that, as is usually the case, 
require diffusion to be expressed in terms of mass fraction gra- 
dients, has been described. 
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Figure 1 | Combined temperature diffusion coefficient for a mixture 
of equal parts argon and nitrogen by mole fraction at 1 atm. The 'without 
correction' curve was calculated including only the first term on the 
right-hand side of equation (25), while the 'corrected' curve was calculated 
using the full expression. 

The first point is an important consideration in the consistent 
definition of the combined diffusion coefficients, and the correct 
choice of stoichiometric coefficients becomes critical when diffusion 
driven by external electric fields is taken into account. 

The second point has to be understood and implemented to cor- 
rectly calculate the combined temperature and pressure diffusion 
coefficients. In particular, very large errors are introduced in the 
calculation of the combined temperature diffusion coefficient if this 
point is neglected. 

Finally, the third point provides an approach to implement the 
combined diffusion coefficients into standard computational models 
of thermal plasmas. 
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